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We revisit the phenomenon of quantum stochastic reso- 
nance in the regime of validity of the Bloch equations. We 
find that a stochastic resonance behavior in the steady-state 
response of the system is present whenever the noise-induced 
relaxation dynamics can be characterized via a single relax- 
ation time scale. The picture is validated by a simple nuclear 
magnetic resonance experiment in water. 
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The interplay between dissipation and coherent driv- 
ing in the presence of dynamical nonlinearities gives rise 
to a variety of intriguing behaviors. The most paradig- 
matic and counterintuitive example is the phenomenon 
of stochastic resonance (SR), whereby the response of 
the system to the driving input signal attains a maxi- 
mum at an optimum noise level Q|. By now, stochastic 
resonance has been demonstrated in several overdamped 
bistable systems as diverse as lasers, semiconductor de- 
vices, SQUID's, and sensory neurons, the required noise 
tuning being accomplished by either controlling the in- 
jection of external noise or by suitably varying the tem- 
perature of the noise- inducing environment. 

Due to the broad typology of situations it can exem- 
plify and its inherent simplicity, a preferred candidate for 
theoretical analysis is represented by a driven dissipative 
two- level system (TLS). The investigation has only re- 
cently been taken into the quantum world, where some 
prominent results have been established for the so-called 
spin-boson model. The latter schematizes the archety- 
pal situation of a driven quantum-mechanical tunneling 
system in contact with a harmonic heat bath, the result- 
ing dissipation being commonly addressed in the linear 
Ohmic regime &-^|- Within this framework, a quantum 
SR phenomenon induced by a resonant irradiation with a 
continuous- wave field has been characterized analytically 
and verified through exact numerical path-integral 
calculations 

In the present work, we show that a stochastic res- 
onance phenomenon occurs for a much wider class of 
driven two-state quantum systems, whose relaxation dy- 
namics can be accounted for by conventional Bloch equa- 
tions. We find that, irrespective of the details of the mi- 
croscopic picture, the essential requirement is the emer- 
gence of a single relaxation time scale. The prediction 
is neatly demonstrated by a nuclear magnetic resonance 
experiment on a water sample. 



Let us consider a two-state quantum system whose 
density operator p is represented in terms of the Bloch 
vector s as p = (l + s-a)/2 i.e., s,(t) = (cr,-(t)), i = 1,2, 3, 
in the customary pseudo-spin formalism J/j. Within the 
semigroup approach for open quantum systems ||, the 
most general (completely) positive relaxation dynamics 
induced by the coupling to some environment is described 
by a quantum Markov master equation of the form 
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-T[H(t),p} + - a k i{[cr k p,(Ti} + [a k ,pai]} 



(1) 
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The Hamiltonian H(t), which describes the interac- 
tion of the TLS with the (classical) driving field, can 
be expressed as H(t) = TiujQO-3/2 + V(t), V(t) = 
h(2u>i) cos(fit)<Ji/2, where the Larmor frequency u>q — 
(E2 — Ei)/h associated with the TLS energy splitting 
and the Rabi frequency ui\ proportional to the alternating 
field amplitude have been introduced. The above Hamil- 
tonian is identical to the one describing a driven tunnel- 
ing process in a symmetric double-well system with "lo- 
calized" states provided by ci-eigenstates: By rotating 
the spin coordinate by 7t/2 about the y-axis, one formally 
recovers the picture of tunneling in the z-rcpresentation 
that is encountered in the literature |]-|L[[(3. The dissi- 
pative component of the TLS dynamics is fully charac- 
terized by the positive-definite 3x3 relaxation matrix 
A = {dki}, determining the equilibrium state of the sys- 
tem and the relaxation time scales connected with the 
equilibration process. The Bloch equations correspond 
to an especially simple realization of A, the non-zero el- 
ements being specified in terms of 3 independent param- 
eters: an - a 22 = (2T 1 )-\ a 33 = (T 2 )- x - (2T!)" 1 , 
ax2 = Q-21 = i(V2Ti)~ 1 s eq . T\,Ti, and s eq are identified 
as the longitudinal and transverse lifetimes, and the equi- 
librium value of the population difference respectively. 
Thus, Eq. (]]]) takes the following familiar form ||: 



s'i = oj s 2 - T 2 1 s 1 
s 2 = -w si - T 2 -1 s 2 
S3 = — 2u)\ cos(f2i)s2 



- 2u>i cos(f2i)s3 , 

rf 1 ^ - Seq) ■ 



(2) 



In microscopic derivations of (|2j), including the ones 
based on the spin-boson model in the appropriate limit 
pO|,p|, relaxation is caused by elementary processes in- 
volving noise-assisted transitions between the TLS en- 
ergy levels or purely dephasing events with no energy ex- 
change between the system and the environment. The 
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overall relaxation rates are obtained by integrat- 
ing such fluctuation and dissipation effects over the en- 
vironmental modes, weighted by the appropriate noise 
spectral densities. Since the latter contain the coupling 
strength between the system and the environment, relax- 
ation rates are themselves directly proportional to the 
underlying noise intensity. Note that the above treat- 
ment in terms of a constant matrix A is only valid for 
external fields that are relatively weak on the TLS energy 
scale i.e., 2\u>%\ -C u>q. In spite of the many restrictions 
involved, it is remarkable that the Bloch equations (|^) 
are of such a wide applicability to cover the majority of 
magnetic or optical resonance experiments. 

For times long compared to the time scales T\ t i of the 
transient dynamics, the motion of the system reaches 
a steady-state behavior that is insensitive to the ini- 
tial condition and acquires the periodicity of the driv- 
ing. In particular, the asymptotic TLS coherence prop- 
erties are captured by the off-diagonal matrix element 
s i(t)i S2(t). It is standard practice to formulate an in- 
put/output problem, where the TLS is regarded as a dy- 
namical system generating si(t) = (ax(t)) as the output 
signal in response to a given input drive V(t). By letting 
limt-»oo (ci (t)) = sf(t) denote the limiting steady-state 
value of si(t), a figure of merit for the system response 
is the so-called fundamental spectral amplitude (!]] , 

f/(n, Wl ) - K(t)\ = &(2wi)|x(fi,wi)| . ( 3 ) 

where the connection to the complex susceptibility 
x(f2, wi) = — ix"(^l,(jJi) is made explicit. 

The competition between driving and dissipative forces 
sets the boundary between the linear vs. nonlinear re- 
sponse regimes. In the limit where lo\ <C T^\, only first- 
order contributions in u)\ are significant and the suscep- 
tibility x(f2, u)x) = x(^) m (H) can be calculated within 
ordinary linear response theory. The linear regime im- 
plies that absorption of energy from the applied field 
occurs without disturbing populations from their equi- 
librium value s eq . Linear behavior breaks down when- 
ever cji > T-j~ 2 . Strongly nonlinear-response regimes can 
be entered for arbitrarily weak fields as long as the cou- 
pling to the environment and the induced noise effects are 
weak enough. For both linear and nonlinear driving, the 
amplitude n(Q, lu\) of the output signal also depends on 
the various parameters characterizing the noise process. 
Quite generally, the phenomenon of stochastic resonance 
can be associated with the appearance of non-monotonic 
dependencies upon noise parameters, leading to the op- 
timization of the response at a finite noise level. 

We focus on the Bloch equations (^) with resonant 
driving, Q = u>q. It is then legitimate to invoke 
the rotating-wave approximation and replace the al- 
ternating field V(t) with V(t) = Tiwi cos(i7t)<7i/2 — 
hu>ism{Vit)a2/2. The rotating-frame description of the 
Bloch vector jl is introduced via the time-dependent rota- 
tion R = exp(iil<j3t/2) i.e., pr = RpR^ 1 = (l+/x-c?)/2, 



p, = (u,v,w) 0. The steady-state solution to the Bloch 
equations is well known jLl 12|. In particular, \' an d x" 
are read from the dispersive and absorptive components 
u, v of the Bloch vector respectively, and the spectral am- 
plitude rj(Q — ujq, uji) = (u 2 +v 2 ) 1 / 2 . A simple expression 
is found for the nonlinear response: 

T](n = U) ,W X ) = S eq .^itrr ■ ( 4 ) 
1 + LU 1 lll2 

Suppose now that we have the capability of manipulat- 
ing the strength of the coupling of the TLS to its environ- 
ment, thereby changing the relaxation times T\, T%, For a 
fixed driving amplitude uj\, rj displays purely monotonic 
behaviors if T\,T% are varied independently. However, 
if a single relaxation time is present, T± = T2 = T12, T] 
develops a local maximum characterized by 



(5) 



The occurrence of such a peak in the steady-state re- 
sponse as a function of the noise strength can be pictured 
as a stochastic resonance effect in the TLS. Physically, 
the condition for the maximum in (|^) can be thought 
of as a synchronization between the periodicity of the 
rotating-frame vector in the absence of relaxation and 
the additional time scale emerging when dissipation is 
present. In semiclassical terms, a constraint of the form 
T\ = kT%, k = const., indicates that the spectral densi- 
ties of the fluctuating environmental fields along different 
directions are not independent upon each other. Simple 
examples include noise processes that effectively originate 
from a single direction or that equally affect the system 
in the three directions. In fact, the existence of a single 
relaxation time is a feature shared with earlier investi- 
gations of quantum SR based on the driven spin-boson 
model J^,^|, where it arises as a necessary consequence 
of the initial assumption that environmental forces ex- 
clusively act along the tunneling axis. However, we em- 
phasize that our discussion is done without reference to a 
specific model, encompassing in principle a larger variety 
of physical situations. 

Apart from this conceptual difference, the SR phe- 
nomenon evidenced above is characterized by the same 
distinctive features found for the spin-boson model on 
resonance |^,^|. According to (^J), the maximum steady- 
state response is independent of the driving amplitude, 
whereas the position of the peak shifts toward shorter 
relaxation times with increasing u)\. Thus, weaker noise 
strengths require weaker input fields to attain a large 
response. This brings the nonlinear nature of the SR 
mechanism to light, for weaker dissipation more eas- 
ily pushes the system into a regime where w\T\2 ^ L 
No SR peak occurs in the limit uj\T\2 <C 1 where lin- 



ear response theory applies and rj 



,(wiTi 2 ) < s eq - 



Thus, SR results in efficient noise-assisted signal ampli- 
fication. Breakdown of linear behavior is more convinc- 
ingly demonstrated by looking at the dependence of the 
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response (^) upon the external field strength. It is easily 
checked that the condition (|J) simultaneously optimizes 
the response against lji, with oj\ = (T1T2) ~ 1//2 . However, 
it is only when T\ = kT2 that the existence of such an 
optimal held amplitude coexists with a SR effect. 

Nuclear Magnetic Resonance (NMR) provides a natu- 
ral candidate for a direct experimental verification of the 
predicted phenomenon. In NMR, the Bloch equations 
describe the motion of the magnetization vector M 
of spin 1/2 nuclei ( 1 H) that are subjected to a static 
magnetic field Bq along the £-axis and a radio-frequency 
signal with amplitude 2\Bi\ <C Bo applied at frequency 
f2 along the i-axis. The mapping is established by iden- 
tifying s = M, s eq = M , uj = jB , uji = 7S1, M and 
7 denoting the equilibrium magnetization and the gyro- 
magnetic ratio respectively. Relaxation processes arise 
due to a multiplicity of microscopic mechanisms [fl2"|| . 
For a liquid spin 1/2 sample, the leading contribution 
arises from fluctuations of the local dipolar field caused 
by bodily motion of the nuclei. The longitudinal relax- 
ation time T\ is essentially determined by the x- and y- 
components of the local magnetic fields at the Larmor 
frequency, while the transverse lifetime T2 takes extra 
contributions from static components of the z-field, im- 
plying that T2 < 2T\ ordinarily |lj,|l^]. Let us assume 
as above that = c^o- Once the steady state is reached, 
the magnetization vector M precesses about the z-axis 
with the periodicity of the r.f. field. The variation of 
the dipole moment in the tranverse plane induces a mea- 
surable e.m.f. in a Faraday coil. This provides access 
to the relevant quantity 77 of Eqs. (|j)-(|]), which repre- 
sents the length of the transverse magnetization vector, 
iu 2 +v 2y/2 = iM 2 {t) + M ^ t) y/2 

Our experiment consists in probing the steady-state 
magnetization response of water as a function of the noise 
strength inducing the natural relaxation processes. The 
1 H Larmor frequency at B = 9.4 T is w /27r = 400 MHz, 
with relaxation times T\ = 3.6 s, T2 = 2.5 s. The sample 
can be brought to a regime where T\ ~ T2 upon addition 
of the paramagnetic salt copper sulfate (CUSO4). With 
concentrations in the range between 40 mM and 100 mM, 
collision events with the impurity dominate the nuclear 
relaxation dynamics. This effectively pushes the system 
into a regime of rapid motion where the correlation time 
of the local magnetic fields seen by the nuclei is very short 
on the scale u/q 1 , thereby ensuring that T\ = T2 — T%% 
. Higher concentrations of the CUSO4 additive result 
in a shorter relaxation time T\% hence implying an effec- 
tive tuning of the noise strength. All measurements were 
performed at room temperature with a Bruker AMX400 
spectrometer on five water samples with additive concen- 
tration in the above range. 

Independent measurements of 2\ and T2 were made to 
confirm that the amount of CUSO4 was sufficient to make 
them equal. T\ was measured via an inversion recovery 
technique O] , by looking at the recovery curve of M z (t) 



after the application of a ir pulse causing M z (0) — —Mq. 
Values of T2 were inferred from the decay of the echo 
signals in a standard Carr-Purcell sequence where it ro- 
tations were used to refocus dephasing due to inhomoge- 
neous broadening JL3| . For the 5 concentrations utilized, 
Ti and T2 were found to be within 1 ms of the average 
value T\2 which is listed for each sample in Table I. 



CuS0 4 (mM) 


T12 (ms) ± 1 ms 


40 


45.5 


50 


36.5 


60 


28.5 


75 


25.0 


100 


18.0 



TABLE I. Relaxation time T12 = T\ = T2 as a function of 
the CuSCU paramagnetic impurity concentration for the 5 
water samples used in the experiment. 

For each sample, the response to a long external r.f. 
pulse was measured for various values of the driving am- 
plitude. For a given driving amplitude, the duration of 
the pulse was increased up to about 200 ms and the 
reading was continued until a constant e.m.f. value was 
reached, confirming that all transient responses had suf- 
ficiently decayed. Under these conditions, the observed 
steady-state value is equivalent to the one produced by 
a cw-irradiation as assumed in Eqs. (0). A delay long 
with respect to T12 was waited between each pulse to 
allow the sample to return to equilibrium. The value of 
the r.f. amplitude was determined by extrapolating mea- 
surements of the nutation rate v\ = lui/2it at a high field 
setting down to the relevant lower-field domain in the 
neighbourhood of ui\ sa T^f . While the relative error be- 
tween two r.f. setting is found to be small, the systematic 
error associated with the extrapolation turns out to be 
significant. A linear correction of the frequency scale was 
included in the analysis to compensate for such error. 

The experimental results are shown in Figs. 1 and 2. 
Fig. 1 evidences the SR peak for three values of the driv- 
ing amplitude. A bell-shaped maximum in the response 
profile is clearly visible, as well as the expected shifting 
of the peak location with increasing uj\. For each curve, 
the SR condition T* 2 « 1 of Eq. (|) is in fairly good 
agreement with the observed behavior, existing discrep- 
ancies being accounted for by the residual error affecting 
the determination of uj±. In Fig. 2 the complementary 
characterization of the SR effect in terms of nonlinear 
response to the driving field is displayed for three val- 
ues of Tyi- In each case, ordinary linear response theory 
is valid for small uj\ to the left side of the maximum. 
Thus, SR reveals itself as a signal optimization marking 
the crossover between linear and nonlinear response. 

Beside validating the predictions from the Bloch equa- 
tions (|^), our experimental results also support the con- 
clusions independently reached in earlier theoretical anal- 
yses P,p[. A few remarks are in order concerning the 
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specific case of NMR. First, the present experiment is 



not a stochastic NMR experiment 14 1. While the ob- 
vious similarity is that both methods probe the nuclear 
spin system by looking at the transverse magnetization 
response, in stochastic NMR the system is directly ex- 
cited by noise, which is therefore always extrinsic (and 
classical) in origin. More importantly, as mentioned al- 
ready, the solutions to the Bloch equation have a long 
history as a tool to investigate magnetic resonance be- 
haviors. In particular, the existence of an optimum r.f. 
amplitude uj"{ — {TiT 2 )~ 1 / 2 is a feature pointed out long 
ago by Bloch himself jj| . However, only the SR paradigm 
provides the motivation to regard relaxation features as 
controllable output parameters and to look at the usual 
response behavior along different axes in the parameter 
space. Even once this is done, this does not automatically 
lead to SR. Rather, it is the recognition that a single axis 
T\ = T2 is effectively needed to bring out the fingerprint 
of the phenomenon and the novel element added to the 
standard NMR analysis. 

In summary, we established both theoretically and ex- 
perimentally the occurrence of stochastic resonance in 
two-state quantum systems whose relaxation dynamics 
are described by Bloch equations. In addition to sub- 
stantially broadening the existing paradigm for stochastic 
resonance in quantum systems, our results point to the 
possibility of characterizing intrinsic relaxation behavior 
via resonance effects. By offering an optimized way for 
input/output transmission against noise, stochastic res- 
onance carries a great potential for systems configured 
to perform specific signal processing and communication 
tasks [Q. In particular, full exploitation of stochastic 
resonance phenomena could potentially disclose a useful 
scenario for reliable transmission of quantum information 
in the presence of environmental noise and decoherence. 
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FIG. 1. Normalized steady-state response rj/s e q vs. re- 
laxation time T12 for resonant driving i}/2ir = 400 MHz at 
different driving amplitudes: Ui/2tt= 6.3 Hz (circles), 5.5 Hz 
(crosses), 4.8 Hz (squares). Solid lines: theoretical predictions 
from Eq. (4), after systematic correction for the frequency (no 
free parameters). 
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FIG. 2. Normalized steady-state response r//s eq vs. Rabi 
frequency u>i/2-k for resonant driving Q/2n = 400 MHz at 
different relaxation times: T12— 18.0 ms (circles), 28.5 ms 
(squares), 45.5 ms (diamonds). Solid lines as above. 
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